Constraints on linear-negative potentials in quintessence and phantom models from 

recent supernova data 
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We study quintessence and phantom field theory models based on linear-negative potentials of 
the form V{(j)) = s <j). We investigate the predicted redshift dependence of the equation of state 
parameter w[z for a wide range of slopes s in both quintessence and phantom models. We use the 
gold dataset of 157 Snia and place constraints on the allowed range of slopes s. We find s = ± 1.6 
for quintessence and s — ±0.7 ± 1 for phantom models (the range is at the 2a level and the units of 
s are in V^MpHo ~ W'^^eV^ where Mp is the Planck mass). In both cases the best fit is very close 
to s ~ corresponding to a cosmological constant. We also show that specific model independent 
parametrizations of w{z) which allow crossing of the phantom divide line w = ~1 (hereafter PDL) 
provide significantly better fits to the data. Unfortunately such crossings are not allowed in any 
phantom or quintessence single field model minimally coupled to gravity. Mixed models (coupled 
phantom-quintessence fields) can in principle lead to a w{z) crossing the PDL but a preliminary 
investigation indicates that this does not happen for natural initial conditions. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



I. INTRODUCTION 

Recent observations have indicated that the universe 
has entered a phase of accelerating expansion (the scale 
factor obeys a > 0) and that the total amount of clus- 
tered matter in the universe is not sufficient for its small 
average spatial curvature. This converging observational 
evidence comes from a diverse set of cosmological data 
which includes observations of type la supernovae 0, Q , 
large scale redshift surveys Q and measurements of the 
cosmic microwave background (CMB) temperature fluc- 
tuations spectrum jj] . The observed accelerating expan- 
sion and flatness of the universe, requires either a mod- 
ified theory of gravity*^ or, in the context of standard 
general relativity, the existence of a smooth energy com- 
ponent with negative pressure termed 'dark energy'^. 
This component is usually described by an equation of 
state parameter w = ^ (the ratio of the homogeneous 
dark energy pressure p over the energy density p). For 
cosmic acceleration, a value oi w < — i is required as 
indicated by the Friedmann equation 
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Even though the cosmological constant remains a viable 
candidate for dark energy, current observational bounds 
0, 111 on the value of the dark energy equation of state 
parameter w(io) at the present time to (corresponding to 
a redshift z = 0) yield w{z = 0) < 
at best fit. 

The role of dark energy can be played by any physical 
field with positive energy and negative pressure which vi- 
olates the strong energy condition p + 3p > {w > ~^). 
Quintessence scalar fields with positive kinetic term 
(—1 < w < — |) violate the strong energy condition but 
not the dominant energy condition p + p > 0. Their 
energy density scales down with the cosmic expansion 



and so does the cosmic acceleration rate. Phantom 
fields ^3 with negative kinetic term (w < —1) violate 
the strong energy condition, the dominant energy condi- 
tion and maybe physically unstable. However, they are 
also consistent with current cosmological data and ac- 
cording to recent studies 0, |1| they maybe favored over 
their quintessence counterparts. 

Homogeneous quintessence or phantom scalar fields are 
described by Lagrangians of the form 



(1.2) 



where the upper (lower) sign corresponds to a 
quintessence (phantom) field in equation (|1.2() and in 
what follows. The corresponding equation of state pa- 
rameter is 
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For quintessence (phantom) models with V{4i) > 
(V{(j)) < 0) the parameter w remains in the range 
— 1 < w < 1. For an arbitrary sign of V{(t>) the above 
restriction does not apply but it is still impossible for w 
to cross the phantom divide line (hereafter PDL) w = — 1 
in a continous manner. The reason is that for w = — 1 
a zero kinetic term zLcj)'^ is required and the continous 
transition from w < —I to w > —1 (or vice versa) would 
require a change of sign of the kinetic term. The sign of 
this term however is fixed in both quintessence and phan- 
tom models. This difficulty in crossing the PDL w = —1 
could play an important role in identifying the correct 
model for dark energy in view of the fact that data fa- 
vor w ~ — 1 and furthermore parametrizations of w{z) 
where the PDL is crossed appear to be favored over the 
cosmological constant w = —1 (see section HI and Refs 

US). 
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In view of the above described problem it is interesting 
to consider the available quintessence and phantom scalar 
field models and compare the consistency with data of 
the predicted forms of w{z) among themselves and with 
arbitrary parametrizations of w{z) that cross the PDL. 
This is the main goal of the present study. 

We focus on a particular class of scalar field potentials 
of the form 



(1.4) 



where we have followed Ref. and set = at F = 0. 
As discussed in section II (see also Ref. the field 

may be assumed to be frozen (0 = 0) at early times due 
to the large cosmic friction H{t). It has been argued ^3 
that such a potential is favored by anthropic principle 
considerations because galaxy formation is possible only 
in regions where V{<j)) is in a narrow range around V = Q 
and in such a range any potential is well approximated 
by a linear function. In addition such a potential can 
provide a potential solution to the cosmic coincidence 
problem [13. 

For quintessence models the scalar field behavior has 
been studied extensively 11, 14, 15, 16J and shown to 
lead to a future collapse of the scale factor (termed 'cos- 
mic doomsday') due to the eventual evolution of the 
scalar field towards negative values of the potential where 
the gravity of the field is attractive. Such a doomsday 
however does not occur in the corresponding phantom 
models because the scalar field moves towards higher val- 
ues of the potential where the field gravity is repulsive 
and leads to faster acceleration of the scale factor (see 
Figure 2 below) and eventually to a Big Rip singular- 
ity Thus w evolves towards values less than —1 for 
phantom models. One of the goals of this paper is to 
compare the consistency with Snia data of this phantom 
behavior of ■u;(z) with the corresponding consistency of 
the quintessence behavior where w evolves towards values 
larger than —1. 

The structure of the paper is the following: In the 
next section we solve numerically the field equations for 
phantom and quintessence models coupled to the Fried- 
man equation and derive the cosmological evolution of 
the scalar field, the scale factor and the equation of state 
parameter w for several values of the potential slope s. 
In section III we fit the derived Hubble parameter to the 
SnIa Gold dataset and obtain constraints for the po- 
tential slope s for both phantom and quintessence mod- 
els. The quality of fit of these models is also compared 
to the quality of fit of arbitrary w{z) parametrizations 
that can cross the PDL. Finally in section IV we sum- 
marize our results and state the main questions that 
emerge from them. In a preliminary effort to address 
some of these issues we show the evolution of w{z) in 
a mixed quintessence -I- phantom model where the dark 
energy consists of a mixture of interacting phantom and 
quintessence fields. 



II. PHANTOM AND QUINTESSENCE FIELD 
DYNAMICS 



In order to study in some detail the scalar field dynam- 
ics, we consider the coupled Friedman-Robertson- Walker 
(FRW) and the scalar field equation 
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where Mp = {8>nG) is the Planck mass and we have 
assumed a potential of the form 
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where the upper (lower) sign corresponds to quintessence 
(phantom) models. By setting 
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ViMpHl ^ ^ 
equation 12.1|l may be written in rescaled form as 



(2.4) 
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while the scalar field equation H2.2|l remains unchanged. 
It is now straightforward to solve the system numerically 
[Tsl using the following initial conditions (t t-i ~ 0) 



a(u) = (^)V3 tr 

HU) = 



(2.6) 



since the universe is matter dominated at early times and 
an infiationary phase would redshift the gradient and ve- 
locity of the scalar field while the large cosmic friction ^ 
would freeze it at early times after inflation. The value 
of (j>i is chosen for each value of the slope s such that 
rio</> — i(p^{to) + V{(j){to j) = 1 — Qom at the present time 
to (defined by a{to) = H{to) — 1). In what follows we 
have assumed a prior of fiom = 0.3. According to the nu- 
merical solution the scalar field is almost frozen at early 
times (when matter dominates) due to the large cosmic 
friction H{t) ~ At approximately the present time 
when the matter density drops and the field potential 
begins to dominate, the lower friction allows the field to 
move down (up) the potential for quintessence (phantom) 
models (see Figure 1). 

When the potential energy dominates, the universe en- 
ters the present accelerating phase. As the field moves 
down (up), the potential energy becomes negative (more 
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agreement with the corresponding results presented here. 
As discussed in the introduction the PDL is not crossed 
for any value of s. Instead, w{z) evolves towards larger 
(smaller) values of —1 for a quintessence (phantom) 
scalar field. 



Quinlessence s- 



FIG. 2; The scale factor evolution for representative (s = 1) 
quintessence and phantom models with linear potential. The 
present time corresponds to to = 0.96 as in Figure 1. 



positive), the field gravity becomes attractive (more re- 
pulsive) and the scale factor begins to decelerate again 
(accelerate more rapidly) until the universe ends with a 
Big Crunch (Big Rip) (see Figure 2). Using the numerical 
solution of the system H2.5|l - H2.2|) we can also evaluate the 
redshift dependence of the equation of state parameter 



w{z) 



p _ ±\^^ - v{4>) _ + 
P ±i(/.2 + y(0) 



(2.7) 



which is shown in Figure 3 in the redshift range < 
z <2 for both quintessence and phantom models and for 
several values of the slope s. 

In the case of quintessence ?£(z) has been evaluated 
in Ref. ^ (see also 0, 0, 0) with results in good 



III. FIT TO THE GOLD DATASET 

Having solved numerically the rescaled system 
(|2.5|I . H2.2|I for both quintessence and phantom models, 
it is straightforward to obtain the corresponding Hub- 
ble parameter H{z] s) — ^{z\ s) as a function of redshift. 
This may now be used to obtain the corresponding Hub- 
ble free luminosity distance 



Dt{z-s) = {l + z) 



dz'^^ 
Hiz';s) 



(3.1) 



Using the maximum likelihood technique |20j we can find 
the goodness of fit to the corresponding observed {zi) 
(i = 1, 157) coming from the Snia of the Gold dataset 
131 . The observational data of the gold dataset are pre- 
sented as the apparent magnitudes m{z) of the SnIa with 
the corresponding redshifts z and la errors <7m(z)- The 
apparent magnitude is connected to Dl{z) as 



m(z; s) = Af(M, Hq) + Uogw{DL{z-, s)) 



(3.2) 



where M is the magnitude zero point offset and depends 
on the absolute magnitude M and on the present Hubble 
parameter Hq as 

-1 



M^M + Uogw{ 



Mpc 



■) + 25 



(3.3) 



The goodness of fit corresponding to any slope s is de- 
termined by the probability distribution of s i. e. 



(3.4) 
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where 

X-(M,.) = f(-°"^--)-""(--^^'-»^ (3.5) 

and A/" is a normalization factor. The parameter M is a 
nuisance parameter and can be marginalized (integrated 
out) leading to a new defined as 

/+00 
e^^ /^dM (3.6) 
-C30 

Using equations l|3.5)l and (|3.6() it is straightforward to 
show (see Refs 8, 21]) that 

X^(s) = x'{M - 0, s) - ^ + ln(C/2^) (3.7) 
where 

Bis) = g (3.8) 

1=1 m°''»(zi) 
157 

C = Y.-2 (3-9) 

i=l 7n"i>^{zi) 

Equivalent to marginalization is the minimization with 
respect to M. It is straightforward to show [23| that 
can be expanded in M as 

X^{s) = x'^iM = 0, s) - 2MB + M^C (3.10) 

which has a minimum for M = ^ at 

x\s)^xHM^O,s)-^ (3.11) 

Using (|3.11|l we can find the best fit value of s (s = Sq) 
as the value that minimizes x^{s) (x^(so) = yl.v.)- The 
1(7 error on s is determined by the relation |2flj 

AxL =x'(si.)-xL„ = 1 (3.12) 

i.e. s is in the range [so,sio-] with 68% probability. 
Similarly the 2a error (95.4% range) is determined by 
AxL = 4 and the 3cr error (99% range) by AxL = 6-63. 

Figures 4 and 5 show plots of the differences Ax^ (s) = 
X^(s) — x^('S — 0) with respect to the cosmological con- 
stant (x^(s ~ 0) = 177.1) for quintessence and phantom 
models with the la, 2a and 3a ranges marked by dashed 
lines. 

From Figures 4 and 5 it is clear that for both phan- 
tom and quintessence models the best fit is obtained for 
s ~ corresponding to a cosmological constant. For 
phantom models however the fit in the range < s < 1 
is almost degenerate (all values of s in this range give es- 
sentially equally good fit) . The 2a range for quintessence 
is s ~ ± 1.6 while for phantom fields the correspond- 
ing range is s ~ ±0.7 ± 1 (the best fit is at s ~ 0.7 and 
the symmetry of the evolution with respect to +s — > —s 
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FIG. 4: The differences Ax^(s) = X^(s) - X^(s - 0) for 
quintessence models. The curve has been sampled at s = 
0, 0.1, 0.2, 0.5, 0.75, 1.0, 1.2, 1.5, 2.0, 2.5, 3.0 and the 
corresponding points have been joined. 
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FIG. 5: The differences Ax^ = X^(s) - X^(s - 0) for 

phantom models. The curve has been sampled at s = 

0, 0.2, 0.5, 0.7, 1.0, 1.5, 2.0, 2.35 and the corresponding 
points have been joined. 



has been imposed). The value of x^ at best fit is 
Xmin — X^(* ~ 0) = 177.1 for quintessence (identical to 
the value corresponding to a cosmological constant with 
flA = 0.7) and xliin = X^(s - 0-7) = 176.9 for phantom 
fields. Clearly both classes of models can not provide 
better fits than the cosmological constant KCDM . 

We are thus faced with the following question: 'What 
are the particular features required by w{z) for better 
fits to the Snia data?'. To address this question we can 
use arbitrary parametrizations of wiz) and identify the 
forms of w{z) that best fit the data. This task has been 
undertaken by several authors '7,'^] and the best fit forms 
of w{z) found had the following common properties: 
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• The value of w{z = 0) at best fit was found to be 
in the range — 1 > ■w{z = 0) > —2 

• The function 'w{z) at best fit was found to cross the 
PDL from below at least once with ^ > in the 
range < z < 1. 

This is demonstrated in Figure 6 where we plot w{z) 
for two representatives of the field theory models studied 
here (quintessence with s — 2 and phantom with s — 1.5) 
superimposed with wlz) for the best fits of two arbitrary 
parametrizations. The paramctrizations considered are 
the following 



A linear ansatz 



— Wo + wi z 



Using the equation^ 



w{z) 



(1 + ^) 



dlnH 



Pde{z) 

pde{z) ~ I - [H^yn^^^ii + z) 



(3.13) 



(3.14) 



we can obtain the Hubble parameter H[z) corre- 
sponding to the w{z) of (|3.13|l as 



H\z)=Hl[nom{l + zf + 

(1 - 17o™)(l + z)3(i+"'«-'^i) e^'^'"-] 



(3.15) 



This can now be used to obtain D'j^{z; wo,wi) from 
equation (|3.1|l and minimize the obtained from 
the Gold datasetQ. Using the Gold dataset, the 
best fit parameter values for this ansatz arefisj 
(wo,wi) = (-1.4±0.1,l-7±0.4) giving = 174.3 
at the minimum (the errors are at the la level). 



The ansatz 



w{z) = Wq + Wi 



l + z 



(3.16) 



which varies between wq at z = and wo + wi 
at z — > oo with crossover at z = 1 where the two 
values contribute equally. The existence of such a 
crossover has the advantage that observations near 
it apply to a reduced parameter phase space, and 
hence the remaining parameter estimates are more 
sensitive 19J. The Hubble parameter corresponding 
to this ansatz is 




0.25 0.5 0.75 i 1.25 1.5 1.75 

z 

FIG. 6: The evolution of w for two representatives of the 
field theory models studied here superposed with two better 
fits obtained by arbitrary parametrization: a linear ansatz 
w{z) — Wo + wi z and the ansatz 117(2:) — wq + wij^. 



refered above {w{z = 0) < — 1 and cross the PDL w — 
-1). 

For each w{z) we have evaluated Ax^ with respect to 
the cosmological constant Xacdm — 

177.1. The field 

theory models which do not cross the PDL have positive 
Ax^ and are therefore worse fits than KCDM . For ex- 
ample a quintessence field with s = 2 gives Ax^ = 7.8 
while a phantom field with s = 1.5 gives Ax^ — 1.4 (see 
Figure 6). In contrast, the arbitrary parametrizations 
which cross the PDL have negative Ax^ and are there- 
fore better fits than KCDM . In particular for the linear 
ansatz (|3.13|l we find Ax^ — —2.8 while for the smoother 
ansatz of (|3.1()|) we find Ax^ — —3.2. These differences 
mean that the point (wo,wi) = (—1,0) corresponding 
to the cosmological constant from the viewpoint of these 
parametrizations, is worse at more than la from the best 
fits obtained from these parametrizations. The problem 
with such w(z) parametrizations is that it seems to be 
highly non-trivial to approximate their behavior using 
field theory models (even exotic ones). 



IV. CONCLUSION - OUTLOOK 



H^{z)=Hf,[n„r,^{l + zf + 

(1 - VLom){l + z)3(i+"'"+«'i) e3'"i(TT7-i)] (3.17) 

The best fit parameter values for this ansatz areflsj 
{'Wo,wi) = (-1.6 ±0.1, 3.3 ±0.5) giving x^ = 173.9 
at the minimum (the errors are at the la level). 

These parametrizations were chosen for their relative 
simplicity and for leading to fairly good fits to the data 
relatively to other parametrizations (see also Ref. |^). 
As seen in Figure 6 they share both of the properties 



We have shown that phantom and quintessence field 
theory models have serious difficulty to exceed the qual- 
ity of fit of a cosmological constant to the Snia Gold 
dataset. In contrast, arbitrary parametrizations of the 
Hubble parameter and of w{z) that cross the PDL can 
perform significantly better in fitting the SnIa redshift 
data. We are thus faced with the question: 'Which field 
theory models can mimic the behavior of the arbitrary 
parametrizations crossing the PDL?' 

As shown in the introduction this type of behavior 
can not be achieved by either quintessence or phantom 



6 



fields. In principle'23^ such behavior could be achieved by 
field theories involving combinations of quintessence and 
phantom fields (called 'quintom' in recent studies [23|). 
However, preliminary investigations have demonstrated 
that mimicking the w{z) behavior indicated by the best 
fit parametrizations requires significant fine tuning even 
for quintom field theories. This can be demonstrated by 
studying the quintom Lagrangian 

= - - S - S 02 - 9 01 02 (4.1) 



which leads to the rescaled dynamical equations 




0.5 1 U 2 

z 



FIG. 7: The quintom model of the Lagrangian 14.11 with 
natural initial conditions and s = 0.2 can mimic quintessence 
and phantom models by changing the sign of the coupling q. 

Even though we have not explored thoroughly the pa- 
rameter space of this quintom model we have shown 
that the natural initial conditions used for the analysis 
of quintessence and phantom systems presented above 
above do not lead to a crossing of the PDL in the way 
indicated by the best fit parametrizations. As shown in 
Figure 7 (obtained with s = 0.2) we can simply mimic 
either the quintessence or the phantom behavior of w{z) 
by changing the sign of the interaction coupling q. 

Alternative field theory models that can in principle 
lead to a w{z) crossing the PDL with a single scalar 
field can be constructed by considering generalizations 
of the kinetic term of the Lagrangian in the form of k- 
essencejEM |2^. An example of such models includes a 
Lagrangian with a field dependent term multiplying the 
kinetic term in the form 

^ = Ifm^ - Viq>) (4.5) 
It has been shown however!?^ than even in such models 
transitions through w = —1 are physically implausible 
because they are either realized by a discrete set of tra- 
jectories in phase space or they are unstable with respect 
to cosmological perturbations. 

The Mathematica 28j file used for the numerical anal- 
ysis and the production of the figures of the paper can 
be downloaded from 18] or sent by e-mail upon request. 
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